Method and apparatus for windowing in entropy encoding

ABSTRACT

The present invention provides efficient window partitioning algorithms for entropy-encoding. The present invention enhances compression performance of entropy encoding based on the approach of modeling a dataset with the frequencies of its n-grams. The present invention may then employ approximation algorithms to compute good partitions in time O(s*log s) and O(s) respectively, for any data segment S with length s.

This application is a continuation of U.S. patent application Ser. No. 10/894,421 filed Jul. 19, 2004 entitled Method And Apparatus For Windowing In Entropy Encoding (currently allowed), now U.S. Pat. No. 7,296,030 which claims the benefit of U.S. Provisional Application Ser. No. 60/487,992 filed on Jul. 17, 2003, each of which is herein incorporated by reference.

The present invention relates generally to data compression and, more particularly, to a method for efficient window partition identification in entropy encoding to enhance compression performance of entropy encoding based on the idea of modeling a dataset with the frequencies of its n-grams.

BACKGROUND OF THE INVENTION

Compression programs routinely limit the data to be compressed together in segments called windows. The process of doing this is called windowing. String-based compression techniques such as Lempel-Ziv or Burrows-Wheeler often use fixed-size windows suitable for in-core processing. Entropy-encoding techniques such as Huffman or arithmetic compression normally do not require windowing except to bound code lengths or to avoid reading large files multiple times. However, these compressors can benefit from windowing when the statistical models change in different file regions. For example, consider a data file made up from four letters in which two letters appear exclusively in the first half of the file while the other two letters appear exclusively in the second half. If all letters appear with the same frequency, a Huffman compressor would normally encode each letter with two bits. On the other hand, each letter can be encoded with a single bit if each half of the file is treated separately. Adaptive techniques such as adaptive Huffman or splay tree do encode data with shifting models but they often produce inferior codes and incur larger costs in both compression and uncompression times than static Huffman.

Therefore, a need exists for a method for efficient window partition identification in entropy encoding, e.g., with performance much better than O(s³) time.

SUMMARY OF THE INVENTION

In one embodiment, the present invention significantly improves the performance of identifying window partitions in entropy encoding. In particular, the present invention, enhances compression performance of entropy encoding based on the idea of modeling a dataset with the frequencies of its n-grams and employs two approximation algorithms to compute good partitions in time O(s*log s) and O(s) respectively, for any data segment S with length s.

BRIEF DESCRIPTION OF THE DRAWINGS

The teaching of the present invention can be readily understood by considering the following detailed description in conjunction with the accompanying drawings, in which:

FIG. 1 illustrates a flowchart of a partitioning method of the present invention for recursively partitioning a given data segment into smaller segments that can be compressed separately;

FIG. 2 illustrates an example of all possible 2-partitions of a data segment P, having length p, with parts having length at least μ(P) long;

FIG. 3 illustrates a flowchart of a faster variation of partitioning method of the present invention for recursively partitioning a given data segment into smaller segments that can be compressed separately; and

FIG. 4 illustrates the present partitioning method implemented using a general purpose computer or any other hardware equivalents.

To facilitate understanding, identical reference numerals have been used, where possible, to designate identical elements that are common to the figures.

DETAILED DESCRIPTION

The present invention relates to data compression using entropy encoding based on the idea of modeling a dataset with the frequencies of its n-grams.

To better understand the present invention, a description of n-grams and its use are first provided. The present invention uses n-grams to model data. For any data segment S with length s, an n-gram in S is a subsequence of n<s consecutive bytes. Assume an arbitrary but fixed n, the notation S_(i) will denote the n-gram in S starting at location i while S[i] denotes the byte at location i. For example, the string S=abababac is of length 8 and has five 4-grams, of which three are distinct: abab, baba and abac. The 4-grams S₀ and S₂ are the same: abab.

In one embodiment, the present invention repeatedly examine n-gram frequencies of given data segments. Thus, it is beneficial if this step can be executed quickly. For any general data segment S, the notation F_(S) shall be used to denote an associative array of frequencies indexed by the n-grams of S. Suppose that F_(S) was initialized to 0's, the below loop computes all n-gram frequencies: for(i=0; i<=s−n; i+=1) F _(S) [S _(i)]+=1;

This loop runs in time O(s) as long as the cost of indexing the array F_(S) can be bounded by some constant. This can be ensured by implementing F_(S) as a hash table indexed by the distinct n-grams. However, hash table look-up cost is significant and the frequency estimates do not always need to be exact. Thus, F_(S) is chosen to be implemented as a normal array of size A by hashing each n-gram S_(i) to an integer via the below hash function with some preselected constant α: χ(S _(i))=(α^(n−1) S[i]+α ^(n−2) S[i]+ . . . +S[i+n−1]) mod A  (Equ. 1)

The above loop then becomes: for(i=0; i<=s−n; i+=1) F _(S)[χ(S _(i))]+=1;

For nontrivial values of n, the loop can be further optimized by exploiting the linearity of the hash function χ to compute χ(S_(i+1)) from χ(S_(i)) via: χ(S_(i+1))={α(χ(S _(i))−α^(n−1) S[i]}+S[i+n]} mod A  (Equ. 2)

Computing frequencies of 1-grams or single letters is, of course, the basis for Huffman and arithmetic coders. For entropy-encoding compressors, n=1 so A=256 would allow F_(S) to index all possible 8-bit bytes at no loss of accuracy. Henceforth, given a data segment S and a frequency array F_(S), it shall be assumed that F_(S) is indexed by mapping the n-grams via the χfunction as described. Therefore, the notation F_(S)[S_(i)] will mean F_(S)[χ(S_(i))].

Entropy-encoding compressors such as Huffman and arithmetic coders compress data based on modeling the probability of symbols which would be 1-grams or bytes. Although these compressors are sensitive to changes in symbol statistics, they often cannot adapt to evolution in the statistical models. Certain adaptive versions of these algorithms can cope with some model changes but tend to produce less efficient codes and have longer running time. Overall, none of these schemes work well when the statistical models abruptly change. For example, Buchsbaum et al. developed a dynamic programming solution to the problem of grouping columns in large tables to enhance compression. In entropy encoding application with any data segment S with s bytes in length, since each byte is treated as a column, the algorithm can be used to compute an optimal partition in O(s³) time. This is too slow for large datasets with size in the megabytes.

To address this criticality, the present invention provides two methods for computing good partitions using approximations, with significant performance enhancements, in o(s*log s) and O(s) time respectively. This also means that there is good potential gain in finding a good partition of the data into sections with sufficiently different symbol statistics using the present invention, then applying the compressors to each section separately.

Let S be a data segment of length s. A partition π of S with k parts divides S into a sequence of consecutive non-overlapping sub-segments (P₁, P₂, . . . P_(k)) that together exactly cover S. A partition with k parts is referred as a k-partition. For a given compressor Γ and a data segment P, let Γ(P) be the compressed length of P. Then, the compressed length of any partition π of S with k parts is:

$\begin{matrix} {{\Gamma(\pi)} = {\sum\limits_{i = 1}^{k}{\Gamma\left( P_{i} \right)}}} & \left( {{Equ}.\mspace{14mu} 3} \right) \end{matrix}$

FIG. 1 illustrates a flowchart of a partition method 100 of the present invention for recursively partitioning a given data segment into smaller segments that can be compressed separately computation complexity of O(slog s) time. Method 100 starts in step 105 and proceeds to step 110.

In step 110, a data segment P is accepted as input to be partitioned. In step 120, the length of segment P, p, is checked if it is smaller than 2*μ(P). If p is smaller than 2*μ(P), then the method terminates in step 160; otherwise, the method proceeds to step 130.

In step 130, among all 2-partitions of P with both partitions having length at least μ(P), find a π such that Γ(π)<Γ(P). In general, for a segment P with length p with parts having length at least μ(P) in length, there will be a total of (p−2μ(P)+1) 2-partitions. FIG. 2 illustrates the possible combinations of 2-partitions for a segment P with parts having length at least μ(P) in length. Each 2-partition contains a left and a right partition. In the first combination, π₁, the left partition have length μ(P) and the right partition have length p−μ(P). In the second combination, π₂, the left partition have length μ(P)+1 and the right partition have length p−μ(P)−1. In the third combination, π₃, the left partition have length μ(P)+2 and the right partition have length p−μ(P)−2. Following this pattern, in the last combination, π_(p−2μ(P)+1), which is the (p−2μ(P)+1)th combination, the left partition have length p−μ(P) and the right partition have length μ(P). Therefore, in step 130, among all possible 2-partition combinations for a segment P with parts having length at least μ(P) in length, the method calculates

${\Gamma\;(\pi)} = {\sum\limits_{i = 1}^{k}\;{\Gamma\left( P_{i} \right)}}$ based on Equ. 3 and check if the condition Γ(π)<Γ(P) is true for each 2-partition.

In step 140, if a 2-partition that meets the condition of Γ(π)<Γ(P), then the method proceeds to step 150. It is possible that more than one 2-partitions, π, can be found to meet such condition. In that case, the choice of a good π is arbitrary and depends on the applications of the algorithm. In one embodiment of the present invention, the good π used can simply be the first π found among all 2-partitions. If a 2-partition cannot be found to meet the condition of Γ(π)<Γ(P) among all 2-partitions, then the method terminates in step 160.

In step 150, for the chosen 2-partition that meet the condition of Γ(π)<Γ(P), the method recursively and independently uses the left partition and the right partition as inputs to method 100.

In one embodiment of the present invention, μ(P) is chosen to be equal to max(p/5, ε), where p is the length of P and ε=2¹², i.e. 4K bytes. This function is used to determine the minimum length of a partition since small data segments compress poorly. By requiring μ(P) to be a fixed fraction of the length of P, the depth of the recursion of method 100 can be bounded by P(log s). Then, the entire algorithm runs in O(γslog s) time where s is the length of the original data segment S and γ is an estimation of the cost to compute the compressed length function Γ on the parts of the candidate partitions.

For a general compressor, the only way to compute Γ might be to invoke the compressor on the data itself and measure the result. In that case γ might be up to O(s). For entropy-encoding compressors, it is possible to define an estimation function with constant time amortized cost. Consider a data segment P of length p at any recursion level in method 100. Let F_(P) the corresponding array of byte frequencies. Shannon's information theory asserts that the number of bits required to encode a byte i with respect to the data in P is log(p/F_(P)[i]) since F_(P)[i]/p is the empirical probability of i. Let X be an estimate for the length of a table of codes or frequencies that a static Huffman or arithmetic compressor would need to decode data. Then, the compressed length of P, Γ_(e)(P), can be estimated with:

$\begin{matrix} {{\Gamma_{e}(P)} = {{\tau + {\sum\limits_{i = 0}^{255}{{F_{p}\lbrack i\rbrack}{\log\left( \frac{p}{F_{p}\lbrack i\rbrack} \right)}}}} = {\tau + {p\;\log\; p} - {\sum\limits_{i = 0}^{255}{{F_{p}\lbrack i\rbrack}{\log\left( {F_{p}\lbrack i\rbrack} \right)}}}}}} & \left( {{Equ}.\mspace{14mu} 4} \right) \end{matrix}$

In one embodiment of the present invention, τ=5b+(256−b) where b is the number of bytes with non-zero frequency. The factor of 5 was chosen because the Huffman encoder used in one embodiment of the present invention guarantees maximal code length 32. The term 256−b estimates the space needed to encode the bytes not appearing in the data, i.e., having zero code length.

Now, let π1=(P₁, P₂) and π2=(Q₁, Q₂) be two 2-partitions of P such that Q₁ is formed by extending P₁ by one byte on the right. Then Q₂ must have been formed by cutting one byte from the left of P₂. Since only a single byte leaves a part or gets added to it, the frequency arrays F_(P) ₁ and F_(P) ₂ can be updated in constant time to form F_(Q) ₁ and F_(Q) ₂ . As a consequence, Γ_(e)(π₂) can be computed in constant time from Γ_(e)(π₁).

Since all 2-partitions of can be generated by exchanging bytes in a loop starting from the partition (φ, P), where φ is a null data segment, step 130 of method 100 can be implemented so that the total running time of all invocations of the compressed length function Γ_(e) is O(p). Thus, the amortized cost of each Γ_(e) is constant. Further, since each recursion level of method 100 only needs two frequency arrays in the computing loop, the required space for method 100 at all recursion levels is bounded by O(log s). Putting everything together, method 100 can be implemented in O(slog s) time and O(s+logs) space where s is the length of the data segment to be partitioned.

For a slight loss in compression performance, it is possible to eliminate the factor log s from the time complexity of method 100. FIG. 3 illustrates a flowchart of a faster partition method 300 of the present invention for recursively partitioning a given data segment into smaller segments that can be compressed separately with computation complexity of O(s) time. Method 300 starts in step 305 and proceeds to step 310.

In step 310, a data segment P is accepted as input to be partitioned. In step 320, the length of segment P, p, is checked if it is smaller than 2*μ(P). If p is smaller than 2μ(P), then the method terminates in step 380; otherwise, the method proceeds to step 330.

In step 330, all 2-partitions with parts having minimum length of μ(P) are first ordered by the length of their left parts. FIG. 2 illustrates an example of the outcome of such an ordering step. In general, for a segment P with length p with parts having length at least μ(P) in length, there will be a total of (p−2μ(P)+1) 2-partitions. In the first 2-partition, π₁, the left partition have length μ(P) and the right partition have length p−μ(P). In the second 2-partition, π₂, the left partition have length μ(P)+1 and the right partition have length p−μ(P)−1. In the third 2-partition, π₃, the left partition have length μ(P)+2 and the right partition have length p−μ(P)−2. Following this pattern, in the last 2-partition, π_(p−2μ(P)+1), which is the (p−2μ(P)+1)th 2-partition, the left partition have length p−μ(P) and the right partition have length μ(P). Then, step 330 initializes the variables i to 1 and N to p−2μ(P).

In step 340, if i is greater than N, then the method terminates in step 380; otherwise, the method proceeds to step 350. In step 350, if Γ(π_(i))<Γ(P) and Γ(π₁₊₁)>Γ(π_(i)), then the method proceeds to step 370; otherwise, the method proceeds to step 360. In step 360, the method increments i by 1 and the proceeds back to step 340. In step 370, the method recursively apply the right partition of π_(i) as input to method 300.

The basic idea behind method 300 is to consider all 2-partitions of S in order starting from (φ, S), where φ is a null data segment. When a partition is found that improves over the encoding of the entire data segment, it is simply split off from its left part, then used to iterate on the rest. The machinery developed earlier to update frequency arrays can be applied straightforwardly here so that method 300 can be implemented in O(s) time and space where s is the length of the data segment to be partitioned.

FIG. 4 illustrates the present partitioning method implemented using a general purpose computer 400 or any other hardware equivalents. For example, the present partitioning methods and data structures can be represented by one or more software applications (or even a combination of software and hardware, e.g., using application specific integrated circuits (ASIC)), where the software is loaded from a storage medium 406, (e.g., a ROM, a magnetic or optical drive or diskette) and operated by the CPU 402 in the memory 404 of the system. As such, the present partitioning methods and data structures of the present invention can be stored on a computer readable medium, e.g., RAM memory, ROM, magnetic or optical drive or diskette and the like.

While various embodiments have been described above, it should be understood that they have been presented by way of example only, and not limitation. Thus, the breadth and scope of a preferred embodiment should not be limited by any of the above-described exemplary embodiments, but should be defined only in accordance with the following claims and their equivalents. 

1. A method for partitioning a data segment, P, with a length p into smaller segments, P_(i), is compressed separately, the method comprising: storing the data segment, P in a computer readable storage medium; dividing the data segment, P into a plurality of 2-partition pairs π where both partitions within each pair have a predefined length of at least μ(P); estimating a compressed length of each partition within each pair; selecting one of the 2-partition pairs π for partitioning the data segment, P, such that r(π)<r(P), where r(π) represents a compressed length of the one of the 2-partition pairs and Γ(P) represents a compressed length of the data segment, P, wherein at least one of: the storing, the dividing, the estimating or the selecting is performed via a processor; and wherein the method has a computation complexity O(plog p).
 2. The method of claim 1, wherein the compressed length of the data segment is estimated using a function of: ${\Gamma_{e}(P)} = {{\tau + {\sum\limits_{i = 0}^{255}{{F_{p}\lbrack i\rbrack}{\log\left( \frac{p}{F_{p}\lbrack i\rbrack} \right)}}}} = {\tau + {p\;\log\; p} - {\sum\limits_{i = 0}^{255}{{F_{p}\lbrack i\rbrack}{\log\left( {F_{p}\lbrack i\rbrack} \right)}}}}}$ where F_(P)[i] is a corresponding array of frequencies for a byte i and τ is an estimate for a length of a table of codes or frequencies that a static Huffman or arithmetic compressor needs to decode data.
 3. The method of claim 1, where ${\Gamma(\pi)} = {\sum\limits_{i = 1}^{k}{\Gamma\left( P_{i} \right)}}$ and μ(P)=max(p/5, ε) and ε is an arbitrary parameter.
 4. The method of claim 1, further comprising: ordering the 2-partition pairs by a length of the left parts of the 2-partition pairs, for all 2-partition pairs of P with parts having a minimum length of μ(P), where μ(P)=max(p/5, ε) and ε is an arbitrary parameter; and finding a π_(i) such that Γ(π_(i))<Γ(P) and Γ(π_(i+1))>Γ(π_(i)).
 5. The method of claim 4, wherein the finding comprises: when the π_(i) is found, then recursively applying a right partition associated with the π_(i).
 6. A computer-readable storage medium having stored thereon a plurality of instructions, the plurality of instructions including instructions which, when executed by a processor, cause the processor to perform a method for partitioning a data segment, P, with a length p into smaller segments, P_(i), is compressed separately, comprising: storing the data segment, P in a computer readable storage medium; dividing the data segment, P into a plurality of 2-partition pairs π where both partitions within each pair have a predefined length of at least μ(P); estimating a compressed length of each partition within each pair; selecting one of the 2-partition pairs π for partitioning the data segment, P, such that r(π)<r(P), where r(π) represents a compressed length of the one of the 2-partition pairs and Γ(P) represents a compressed length of the data segment, P; and wherein the method has a computation complexity O(plog p).
 7. The computer-readable storage medium of claim 6, wherein the compressed length of the data segment is estimated using a function of: ${\Gamma_{e}(P)} = {{\tau + {\sum\limits_{i = 0}^{255}{{F_{p}\lbrack i\rbrack}{\log\left( \frac{p}{F_{p}\lbrack i\rbrack} \right)}}}} = {\tau + {p\;\log\; p} - {\sum\limits_{i = 0}^{255}{{F_{p}\lbrack i\rbrack}{\log\left( {F_{p}\lbrack i\rbrack} \right)}}}}}$ where F_(P)[i] is a corresponding array of frequencies for a byte i and τ is an estimate for a length of a table of codes or frequencies that a static Huffman or arithmetic compressor needs to decode data.
 8. The computer-readable storage medium of claim 6, where ${\Gamma(\pi)} = {\sum\limits_{i = 1}^{k}{\Gamma\left( P_{i} \right)}}$ and μ(P)=max(p/5, ε) and ε is an arbitrary parameter.
 9. The computer-readable storage medium of claim 6, further comprising: ordering the 2-partition pairs by a length of the left parts of the 2-partition pairs, for all 2-partition pairs of P with parts having minimum length of μ(P), where μ(P)=max(p/5, ε) and ε is an arbitrary parameter; and finding a π_(i) such that Γ(π_(i))<Γ(P) and Γ(π_(i+1))>Γ(π_(i)).
 10. The computer-readable storage medium of claim 9, wherein the finding comprises: when the π_(i) is found, then recursively applying a right partition associated with the π_(i).
 11. An apparatus for performing a method for partitioning a data segment, P, with a length p into smaller segments, P_(i), is compressed separately, comprising: means for storing the data segment, P in a computer readable storage medium; means for dividing the data segment, P into a plurality of 2-partition pairs π where both partitions within each pair have a predefined length of at least μ(P); means for estimating a compressed length of each partition within each pair; means for selecting one of the 2-partition pairs π for partitioning the data segment, P, such that r(π)<r(P), where r(π) represents a compressed length of the one of the 2-partition pairs and Γ(P) represents a compressed length of the data segment, P; and wherein the method has a computation complexity O(plog p).
 12. The apparatus of claim 11, wherein the compressed length of the data segment is estimated using a function of: ${\Gamma_{e}(P)} = {{\tau + {\sum\limits_{i = 0}^{255}\;{{F_{p}\lbrack i\rbrack}{\log\left( \frac{p}{F_{p}\lbrack i\rbrack} \right)}}}} = {\tau + {p\;\log\; p} - {\sum\limits_{i = 0}^{255}\;{{F_{p}\lbrack i\rbrack}{\log\left( {F_{p}\lbrack i\rbrack} \right)}}}}}$ where F_(P)[i] is a corresponding array of frequencies for a byte i and τ is an estimate for a length of a table of codes or frequencies that a static Huffman or arithmetic compressor needs to decode data.
 13. The apparatus of claim 11, where ${\Gamma(\pi)} = {\sum\limits_{i = 1}^{k}\;{\Gamma\left( P_{i} \right)}}$ and μ(P)=max(p/5, ε) and ε is an arbitrary parameter.
 14. The apparatus of claim 11, further comprising: means for ordering the 2-partition pairs by a length of the left parts of the 2-partition pairs, for all 2-partition pairs of P with parts having a minimum length of μ(P), where μP)=max(p/5, ε) and ε is an arbitrary parameter; and means for finding a π_(i) such that Γ(π_(i))<Γ(P) and Γ(π_(i +1))>Γ(π_(i)).
 15. The apparatus of claim 14, wherein the finding means, when the π_(i) is found, then recursively applies a right partition associated with the π_(i). 